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The figure of merit for refrigerators performing finite-time Carnot-like cycles between two reser- 
voirs at temperature Th and T c (< Th) is optimized. It is found that the coefficient of performance 
at maximum figure of merit is bounded between and (-y/9 + 8e c — 3)/2 for the low-dissipation 
refrigerators, where e c = T c / (Th — T c ) is the Carnot coefficient of performance for reversible re- 
frigerators. These bounds can be reached for extremely asymmetric low-dissipation cases when the 
ratio between the dissipation constants of the processes in contact with the cold and hot reservoirs 
approaches to zero or infinity, respectively. The observed coefficients of performance for real refrig- 
erators are located in the region between the lower and upper bounds, which is in good agreement 
with our theoretical estimation. 

PACS numbers: 05.70.Ln 



Introduction.- The issue of efficiency at maximum 
power output has attracted much attention since the 
seminal achievements made by Yvon [l|, Novikov 
Chambadal "S\ , Curzon and Ahlborn 4] , which gives rise 
to finite-time thermodynamics, a new branch of non- 
equilibrium thermodynamics, and opens open new av- 
enues to the perspective of establishing more realistic 
theoretical bounds for real heat engines as well as re- 
frigerators d-Q. 

Previous reported works on this subject show that dif- 
ferent model systems exhibit various kinds of behaviors 
at large relative temperature difference between two ther- 
mal reservoirs at temperatures Th and T c (< Th), in spite 
that they show certain universal behavior at small rela- 
tive temperature difference f9j-fl7j leading to recent dis- 
cussions on the bounds of efficiency at maximum power 
output for Carnot-like heat engines (l8l - [23j . In particu- 
lar, Esposito et al. investigated low-dissipation Carnot- 
like engines by assuming that the irreversible entropy 
production in each isothermal process is inversely pro- 
portional to the time required for completing that pro- 
cess [l||. Furthermore, they obtained that the efficiency 
at maximum power output for low-dissipation engines is 
bounded between ?y_ = r/c/2 and 7]+ = r)c/(2 — T)c) P^t - 
where r\c = 1 — T c /Th is the Carnot efficiency of re- 
versible heat engines. Besides, Ma [24| proposed the per- 
unit-time efficiency to be another criterion, which can be 
viewed as a compromise between the efficiency and the 
speed of the whole cycle. Two of present authors and 
their coworkers [25[ proved that the efficiency of endore- 
versible heat engines performing at maximum per-unit- 
time efficiency is bounded between rjc/2 and 1 — VI — t]c- 

However, it is relatively difficult to define an optimal 
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criterion and obtain its corresponding coefficient of per- 
formance (COP) for refrigerators [26-35] in the way as 
we address the issue of efficiency at maximum power for 
heat engines provided that minimum power input is not 
an appropriate figure of merit in Carnot-like refrigera- 
tors. Velasco et al. [26| adopted the per-unit-time COP 
as an target function and proved Sca = V £ c + 1 — 1 to 
be the upper bound of COP for endoreversible refrigera- 
tors operating at the maximum per-unit-time COP, being 
e c = T c /(T h - T c ) the Carnot COP for reversible refrig- 
erators. Allahverdyan et al. [13] investigated a quantum 
model which consists of two n-level systems interacting 
via a pulsed external field and took eQ c as the target 
function, where e and Q c are the COP of refrigerators 
and the heat absorbed from the cold reservoir, respec- 
tively. They also proved that the COP of this model 
is bounded between eqa and e c at the small relative 
temperature difference. Chen and Yan [28| suggested to 
take x = zQc/tcycic as the target function, where i cyc i G 
is the time for completing the whole Carnot-like cycle. 
Recently, de Tomas and two of present authors [29| opti- 
mized x for symmetric low-dissipation refrigerators and 
derived the COP at maximum x to be Eca — \/l + £ c— 1- 
The above results give rise to two straightforward ques- 
tions: (i) What target function could be appropriate as 
the figure of merit for refrigerators? (ii) Can we derive 
the bounds of COP at maximum figure of merit for gen- 
eral low-dissipation refrigerators as a counterpart to the 
bounds of efficiency at maximum power output for heat 
engines? We will address these problems in this work. 
We select x = eQc/icyde as the figure of merit and de- 
rive that the COP at maximum figure of merit is bounded 
between and (V9 + 8e c — 3)/2 for low-dissipation refrig- 
erators. Our theoretical prediction is in good agreement 
with the observed data from real refrigerators, which sug- 
gests x — e Qcl '^cycic is appropriate as the figure of merit 
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for refrigerators. 

Model— The refrigerator that we consider performs a 
Carnot-like cycle consisting of two isothermal processes 
and two adiabatic steps as follows. It must be noted that 
the word "isothermal" in this work also merely indicates 
that the working fluid is in contact with a reservoir at 
constant temperature. Here we do not introduce the ef- 
fective temperature of working fluid because the effective 
temperature might not be well-defined in many cases [23[ ■ 

(SI) Isothermal expansion. The working substance is 
in contact with a cold reservoir at temperature T c and 
the constraint on the system is loosened according to 
the external controlled parameter A c (r) during the time 
interval < r < t c , where r is a time variable. It is in the 
sense of loosening the constraint that this step is called 
an expansion process. A certain amount of heat Q c is 
absorbed from the cold reservoir. Then the variation of 
entropy in this process can be expressed as 



AS C = Qc/T c + AS l c 



(1) 



where AS* 1 * > is the irreversible entropy production. 
We adopt the convention that the heat absorbed by the 
refrigerator is positive, so AS** 1 " < AS C . 

(52) Adiabatic compression. This step is idealized as 
the working substance suddenly decouples from the cold 
reservoir and then comes into contact with the hot reser- 
voir instantaneously at time t c . During this transition, 
the controlled parameter is switched from X c (t c ) to Xh(th) 
[> A c (t c )], that is, the constraint on the system is en- 
hanced. It is in the sense of enhancing the constraint 
that this step is called a compression process. There is 
no heat exchange in this transition, i.e. Q2 = 0. The 
distribution function of molecules of working substance 
is unchanged. Thus there is no entropy production in 
this transition, i.e. AS2 = 0. 

(53) Isothermal compression. The working substance 
is in contact with a hot reservoir at temperature Th and 
the constraint on the system is further enhanced accord- 
ing to the external controlled parameter A^ (r) during the 
time interval t c < t < t c + th- A certain amount of heat 
Qh is released to the hot reservoir T/,. Thus the total 
variation of entropy in this process is 



AS h = -Q h /T h + as; 



h - 



(2) 



where AS h r > is the irreversible entropy production. 

(S4) Adiabatic expansion. Similar to the adiabatic 
compression process, the working substance suddenly 
decouples from the hot reservoir and then comes into 
contact with the cold reservoir instantaneously at time 
t c + th- During this transition, the controlled parame- 
ter is switched from Xh(t c + th) to A c (0) [< Xh(t c + th)], 
that is, the constraint on the system is loosened. In this 
transition, both the heat exchange and the entropy pro- 
duction are vanishing, i.e. Q4 = and AS4 = 0. 

Optimizing the figure of merit. -Having undergone a 
whole cycle, the system recovers its initial state. Thus 



the change of entropy is vanishing for the whole cy- 
cle, from which we can easily derive that the varia- 
tions of entropy in two "isothermal" processes satisfy 
AS C = —ASh = AS > 0. Similarly, the total energy also 
remains unchanged for the whole cycle, thus the work in- 
put in the cycle can be expressed as W = Qh — Q c , and 
then the COP of refrigerators is reduced to 



e = Qc/(Qh - Q c ). 



(3) 



Considering Eqs. (P)-© and i cyc i c = t c + th, the figure 
of merit x = ^Qc/t C ycie is transformed into 



X 



T C 2 (A<S — AS, 



ir\2 



[{T h - T C )AS + T C AS? + T h AS h r ](t h + t c ) ' 



(4) 



The variation of entropy AS is a state variable only de- 
pending on the initial and final states of the isother- 
mal processes while AS l c r and AS h r are process vari- 
ables relying on the detailed protocols A(r). In addi- 
tion, ASl r < AS according to Eq. Q. Thus Eq. g]) 
implies that the maximum of the figure of merit cor- 
responds to minimizing irreversible entropy production 
AS l c r and ASjf with respect to the protocols for given 
time intervals t c and th, which is equivalent to that ob- 
tained for Carnot-like heat engines working at maximum 
power output. 

To continue our analysis, we denote the minimum ir- 
reversible entropy production with the optimized proto- 
cols as min{A5f } = L c (t c ) and mm{AS h r } = L h {t h ). 
Intuitively, L c {t c ) and Lh(th) are the monotonous de- 
creasing functions of t c and th, respectively, because the 
larger time for completing the isothermal steps, the closer 
these steps are to quasistatic processes so that the irre- 
versible entropy production ASl r and AS h r become much 
smaller. In particular, AS l c r and AS' l /[ should vanish in 
the long-time limit t c — > 00 and th — > 00. For conve- 
nience, we can make a variable transformation x c = l/t c 
and Xh = l/th- If we consider Eqs. |T]) and ([2]), the heat 
Q c and Qh can be expressed as 



and 



Q C = T C [AS-L c (x c )}, 



Qh = T h [AS + Lh{x h )]. 



(5) 



(6) 



Substituting Eqs. (|5|) and ((6|) into ((3]), we derive the 
COP of refrigerators to be 



T C {AS - L c ) 



Qh - Qc {T h - T C )AS + T C L C + T h L h 



(7) 



Considering i cyc i c = t c + th = l/x c + 1/xh and the 
above equations ©-((7]), we optimize the figure of merit 
X — eQcAcycic with respect to Xh and x c and derive the 
following two equations: 

{Qh - Qc)x h = {IQhIQc - l)T c L' c x c {x h + x c ), (8) 
{Qh - Qc)x c = T h L' h Xh{x h + x c ), (9) 
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where L' c = dL c /dx c and L' h = dL^/dxh- 

Considering Eqs. ©-O and then dividing Eq. flSJ) by 
Eq. (j9|), we can derive that the COP at maximum figure 
of merit satisfies 

e*T h L' h xl = (e*+2)T c L' c x 2 c . (10) 

Similarly, adding Eq. dHJ and Eq. ((9]), we can derive 



1 _ 1 1 

~* ~ ~c + Ne, + (2e c -e,)M/(l + e c ) 



(11) 



with reducing parameters N = (L' c x c + L' h Xh)/(L c + Lh), 
M = L' c x c /(L C + L h ) and e c = T c /(T h - T c ). 




FIG. 1. (Color online) Schematic diagrams of function. The 
dashed line is the diagram of f(y) = 2e c y 2 — 3y — 1 while 
the solid lines correspond to diagram of g(y) = a(l + 2j/) 1 ^ 2 
with different values of a. The points of intersection represent 
the solutions to Eq. I|14[) for different values of a. yo is the 
solution to f(y) = 2e c y — 3y — 1 = 0. 

Bounds of COP at maximum figure of merit-Now we 
turn to the low-dissipation refrigerators by assuming that 
L' c = E c and L' h = E^ are two dissipation constants as 
Esposito et al. [19[ proposed for low-dissipation heat en- 
gines. In this case, N = 1 and M = Y, c x c /(Y< c x c + Y,f l X} l ). 
Particularly, E c = E/, = E for the symmetric low- 
dissipation cases investigated by de Tomas et al. [29j j . 
it is not hard for us to recover its COP at maximum 
maximum figure of merit to be £ca = \A + e c — 1 from 
Eqs. I|10p and Ijlip. However, for the asymmetric low- 
dissipation cases where E c ^ E^, it is more difficult to 
obtain a concise analytic expression of e* than the sym- 
metric case. But we can still estimate its bounds from 
Eq. (jTTj) . According to this equation, we have 



r c h/l + 8(l + e c )/M-3] 
2[(l + e c )/M-l] 



(12) 



which is the key equation in the present work. Because 
e c > and < M < 1, it is easy to prove that e* is a 
monotonous increasing function of M. As a main result, 
from Eq. (|12p we obtain the wished bounds as: 



It is noted that M is also constrained by Eq. (|10l) . 
which pushes us to further discuss the accessibility of 
the lower bound e_ = and the upper bound e+ = 
(V9 + 8e c - 3)/2. Eliminating x c /xh from Eqs. (JTUJ) and 
(fTTj). we have 



2e c y 2 - 3y - 1 = a(l + 2y) 



1/2 



(14) 



where y = 1/e* and a = y / Ti l ^h/T c 'S c . In Fig. [I] we 
schematically plot the function f{y) — 2e c y 2 — 3y — 1 
(dashed line) and g(y) = a(l + 2y) 1 / 2 for different values 
of a (solid lines). The points of intersection between the 
dashed line and solid lines correspond to the solutions to 
Eq. (fT4| for different values of a. It follows that the solu- 
tions to Eq. (|14j) increase with the increasing value of a. 
On the other hand, the values of a can be taken from to 
oo. Therefore we infer that the solutions to Eq. (fl4l) are 
between y = (y/9 + 8e c + 3)/4e c [solution to Eq. flH]) for 
a = 0, i.e. E c /E/i —> oo] and oo [solution to Eq. (fl4")l for 
q4oo, i.e. E c /E/j — > 0]. Noting that y = 1/e*, we ar- 
rive at < < l/yo = (%/9 + 8e c — 3)/2 which is exactly 
the same as inequality (fT3]) . Simultaneously, we obtain 
the condition for reaching the lower and upper bounds: 
e* -> when E c /E h -> and e* -> (-y/9 + 8e c - 3)/2 
when E c /Eft — >• 00. That is, the lower and upper bounds 
of COP at maximum figure of merit can be reached for 
extremely asymmetric low-dissipation refrigerators. Al- 
though the lower and upper bounds of efficiency at max- 
imum power output can also be reached for extremely 
asymmetric low-dissipation heat engines [19j . the subtle 
difference is that the lower bound can be reached when 
E c /E/j — > 00 while the upper one can be reached when 
E c /E/j — > 0, which is in the inverse situation with respect 
to the refrigerators. However, this difference is quite 
reasonable because refrigerators need the input work to 
pump heat from the cold reservoir while heat engines uti- 
lize heat from the hot source to generate work. 




FIG. 2. (Color online) Numerical solutions to Eq. (TTiJ). The 
used values of parameter E c /Eh are marked nearby each 
curves. 



< e* < (V9 + 8e c - 3)/2. 



(13) 



The numerical solutions to Eq. (|T4| can also be calcu- 
lated by setting different values of ratio E c /E^. The cor- 
responding values of e* = 1/y are shown in Fig.(2J from 
which we find that the COP at maximum figure of merit 
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indeed reaches the upper bound e + = (y/9 + 8e c — 3)/2 
when the ratio £ C /Eh is relatively large while it ap- 
proaches the lower bound e_ = when the ratio S c /S/j 
is small enough. In addition, the curve with parameter 
S c /Sfi = 1 corresponds to &ca = \A + £ c — 1, which is 
also located in the region bounded between e_ and £+. 




FIG. 3. (Color online) Comparison between theoretical pre- 
diction and the observed COPs of refrigerators. The circles 
represent the relationship between observed COPs and the 
Carnot COPs calculated according to the working tempera- 
ture region for the reciprocating chiller with nominal cool- 
ing rate 1172 kw while the squares represent that for the 
water cooled reciprocating chiller with nominal cooling rate 
10.5 kw [36|| . The solid line represents the theoretical upper 
bound e+ = (V9 + 8e c - 3)/2. 

Now we compare our prediction with the observed 
COPs of some real refrigerators. The circles and squares 
in Fig. [3] respectively show the relationship between the 
observed COPs of two different kinds of real refriger- 
ators working in different temperature regions and the 



corresponding Carnot COPs calculated according to the 
working temperatures. The raw data are adapted from 
Tables 6.1 and 10.1 in Ref. [Hj]. We stress from this 
figure that all data are located between the optimized 
COPs at E c /S ft -> oo (the solid line) and T, c fS h = 0.01 
(the dotted line), which reveals the capability of the low- 
dissipation assumption and the bounds of the optimized 
COP in order to reasonably estimate the experimental 
results for real refrigerators. Additionally, we also plot 
£ca — \/l + e c — 1 as the dashed line in Fig. [31 from 
which we see that Eqa is neither the upper bound nor 
lower bound of observed COPs. This result suggests that 
X — zQc/tcycic is indeed a very valuable figure of merit 
in comparing with experimental refrigerators data. 

Conclusion - The issue of COP at maximum figure of 
merit for Carnot-like refrigerators is addressed. We ob- 
tain the universal lower and upper bounds of COP at 
maximum figure of merit for low-dissipation Carnot-like 
refrigerators. These bounds can be reached for extremely 
asymmetric dissipation cases. We compare our predic- 
tion with the observed COPs of real refrigerators and 
find that all measured COPs are located in between the 
prediction model. From a theoretical point of view, these 
results for low-dissipation refrigerators can be regarded 
as a counterpart of the bounds of efficiency at maximum 
power output obtained by Esposito et at. [l|| for low- 
dissipation heat engines. In the future work, we will ex- 
tend our discussions to the refrigerators working out of 
the low-dissipation regime based on the key equation (|12p 
and our previous investigation on heat engines 137] . 
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